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We investigate the prediction on the lepton phases in theories with a non trivial correlation 
between quark (CKM) and lepton (PMNS) mixing matrices. We show that the actual evidence, 
under the only assumption that the correlation matrix V product of CKM and PMNS has 
a zero in the entry (1,3), gives us a prediction for the three CP- violating invariants J, Si, and 
S2. A better determination of the lepton mixing angles will give stronger prediction for the CP- 
violating invariants in the lepton sector. These will be tested in the next generation experiments. 
To clarify how our prediction works, we show how a model based on a Grand Unified Theory and 
the permutation flavor symmetry S3 predicts V13 = 0. 
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I. INTRODUCTION 



After the recent experimental evidences about neutrino 
physics PI -[H| we know very well almost all the param- 
eters in the quark [l5[ and lepton [ll|-[35| sectors. We 
measured all the quark and charged lepton masses, and 
the value of the difference between the square of the neu- 



2 



mh and Srni 
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trino masses 5mf 2 
We also know the value of the quark mixing angles and 
phases, and the two mixing angles #12 and 023 in the lep- 
ton sector. The challenger for the next future [3r3|-[53j 
will be to determine the sign of Sm\^ (i.e. the hierarchy 
in the neutrino sector), the absolute scale of the neutrino 
masses, and the value of the 3rd lepton mixing angle #13 
(in particular if is it zero or not). Finally, if 613 is not 
too small, there is a hope to measure the CP violating 
phases. 

Despite the fact that we have all these experimental 
evidences, still a theory of flavor is missing. The main 
point is related to the fact that there is a hierarchy in the 
mixing angles and masses. In particular the quark mixing 
angles are in general much smaller then the correspond- 
ing lepton mixing angles. The small quark mixing angles 
can be explained with continuous flavor symmetries (see 
for instance [54j|). while the structure of the lepton mix- 
ing matrix seems to be better accommodated in models 
with discrete symmetries (sec for example [55]- 58]). Re- 
cently, the disparity that nature indicates between quark 
and lepton mixing angles has been viewed in terms of 
a 'quark-lepton complementarity' (QLC) [H, H(| which 



can be expressed in the relations 



qPMNS 1 nCKM 



45° 



qPMNS 1 nCKM 
'23 > p 23 



45°. (1) 



These relations are related to the parametrization used 
for the CKM and PMNS mixing matrix. From a more 
general point of view, we can define a correlation matrix 
V M as the product of the PMNS [HHH and CKM [H, 
IfJij ] mixing matrices, 



V 



M 



CKM jjPMNS 



(2) 



A lot of efforts have been done to obtain the most fa- 
vorite pattern for the matrix V M [^-[lH. The naive 
QLC relations in eq. |T]) seems to implies V M to be Bi- 
Maximal, i.e. in the standard parametrization it contains 
two maximal mixing angle , and a third angle to be zero. 
A BiMaximal Vm implies a nice, in the sense that it is 
testable at the near future neutrino experiments, predic- 
tion on the 0i^ INS mixing angle [59| and the CP violat- 
ing parameters for the lepton sector [7^]. At first order 
approximation however V M BiMaximal seems not to be 
compatible with the experiments [6 71 ] . From our previous 
work [65[ we learn that V M BiMaximal, although it is not 
ruled out by the experiments, is excluded at 90% CL in 
non SUSY models, or in SUSY models with tan/3 < 40 
where the RGE correction are negligible [73l|-[7r|. 

Despite the fact that the correlation matrix V M cannot 
be BiMaximal, still there are very nice phenomenologi- 
cal consequences from a non trivial V M . All of these 
consequences seems to be related to the fact that exper- 
imental evidences tell us 
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that V$ 
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= is in good 
agreement with the experimental data, and is even the 
preferred value. From a theoretical point of view we like 
very much the fact that experimental data tell us that 
V13 — 0. First of all the theoretical ingredient of the 
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quark-lepton complementarity that gives phenomenolog- 
ical predictions is = 0. In fact V-^ = 0, without any 
other assumptions on the full Vm matrix with the excep- 
tion of the compatibility with the actual experimental 
evidences, implies a testable prediction for the undeter- 
mined lepton mixing angle 

9[ 3 MNS = (9±l)° (see [65]) 

Second because V^ 1 = can be achieved in a natural 
way in theoretical models by imposing some symmetry. 

In this work first we analyze the consequences of 
V\i = on the CP violating invariants for the lepton 
sector. Then we show how a toy model, based on S3 fla- 
vor permutation symmetry (77| and GUT gives V-^ = 0. 

We use a Monte Carlo simulation with two-sided Gaus- 
sian distributions around the mean values of the observ- 
ables to extract the J, Si, and S2 invariants. The in- 
put information on 0f^ INS is taken from the analysis of 
ref. [H| which uses quark lepton correlation (V^ = 0), 
and neutrino and quark data. We obtain that the main 
ingredient to obtain a prediction on the CP violating in- 
variants for the lepton sector is the constraint = 0. 

From the theoretical point of view we show that, if we 
do not take care about a required fine-tuning in obtaining 
the masses, it is possible to construct a toy model [77j . 
based on S3 flavor permutation symmetry and GUT, that 
gives us V$ = 0. In Grand Unified Theories (GUT), 
as long as quarks and leptons are inserted in the same 
representation of the gauge group, we need to include 
in the definition of V M non trivial phases between the 
CKM and PAINS mixing matrices. In our model V M 
is defined by 

t/M _ tjCKM ttPMNS (q\ 



where is a diagonal matrix = diag(e IWi ) and the three 
phases u>i are free parameters (in the sense that they are 
not determined by present experimental evidences). The 
matrix Vm is related to the Dirac and Majorana neutrino 
mass matrix. For this reason its form is given by the 
symmetries of the model. 

The paper is organized as follows: in section [II] we in- 
troduce our notation and the parameterization for CKM 
and PMNS mixing matrices. With the aid of a Monte 
Carlo simulation, we study the numerical correlations of 
the lepton CP violating phases J, Si, and S2 with re- 
spect to the mixing angle &12 . In section IIIII we 
clarify the relation between the correlation matrix V M 
and the neutrino mass matrices in GUT model, with a 
see-saw of type I. After that we show how a S3 flavor 
symmetry can give us V-fif = 0. Finally in section TlVI we 
present a summary and our conclusions. 



II. CP VIOLATING INVARIANTS IN THE 
LEPTON SECTOR 

As usually, we parameterize the lepton mixing matrix 

as 

jjPMNS _ ^23^13^ £/12$m ^ 

where $ is a diagonal matrix with elements {1,1, e*^}, 
and <j> is the Dirac CP violating phase. <I> m contain the 
Majorana phases and is a diagonal matrix with elements 
given by {e 1 ^ 1 , e lc ^ 2 , 1}, and ?7 y are rotation in the 
plan. There are two kind of invariants parameterizing 
CP violating effect. The Jarlskog invariant [28] J that 
parametrizes the effects related to the Dirac phase, and 
the two invariants Si and S2 that parametrize the ef- 
fects related to the Majorana phases. The J invariant 
describes all CP breaking observables in neutrino oscil- 
lations [29] . It is the equivalent of the Jarlskog invariant 
in the quark sector. It is given by 

J = Im{U v ^ U^ 2 U: eU2 U*^ } . (5) 

In the parametrization of eq. (|4]) one has 

J = - sin 2#i2 sin 2^23 sin2#i 3 cos #13 sin</>. (6) 
8 

Then we have the two invariants Si and S2 that are re- 
lated to the Majorana phases. They are 

51 = Im{U v ^K^} (7) 

5 2 = Im{U VeV JJ^ V3 } 

In the parametrization of eq. (|4]) we have 

51 = - cos 0i2 sin 26*13 sin(0 + <p\ ) (8) 

5 2 = - sin #12 sin 2#i 3 sin(0 + <p 2 ) 

The two Majorana phases appear in Si and S2 but not 
in J. 



A. Prediction for CP violating invariants 

In this section we investigate the consequences of a V M 
correlation matrix with a zero (1,3) entry on the unde- 
termined parameters J, Si, and S2. We remember that J 
is the Dirac invariant CP-violating phase, and is the only 
one that can be observed in neutrino oscillations experi- 
ments. As shown in a previous paper [651 ] . the data favors 
a vanishing (1, 3) entry in the correlation matrix V M . So 
in the whole analysis we fix sin 2 9X^' = 0. Moreover 
tan 2 9\" and tan 2 #23 are allowed to vary respectively 
within the intervals [0.3, 1.0] and [0.5, 1.4]. We introduce 
the unitary Wolfenstein parameterization in terms of the 
variables A ,A, p, 77 [78| 

jjCKM = u2 3 $ jjU $ f [/12 5 (9) 
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where one has the relations [79( 

sin^ M = A 
sinfl™ = AX 2 



sin 



■)CKA-r-iS CKM 
'13 e 



AX 3 (p-i v ) 



to all orders in A. We use the updated values for the 
CKM mixing matrix, given at 95%CL by [l5[ 



A = 0.2265t°-221? , A = 0.801 +0 ;ii;; ' 
Tj = 0.189!£iE , p = 0.358_ 0lls , . 



-0.0041 
-,+0.182 



0.041 I 
+0.086 



with 



p + irj = 



VI - A 2 X i (p + irj) 



VT=J?[1- A 2 X i (p + ir])] 
For the lepton mixing angle we impose [HI [l7 
S in 2 ^ 3 MJVS ^ 0.44 X (Itfcg) 



(10) 



(11) 



sin 



2 nPMNS 
7 12 



= 0.314 x (l±8il) , 



and |6E 



qPAINS 
'13 



= 9 



(12) 



(13) 



We allow the JJ CKM parameters to vary, with a two-sided 
Gaussian distribution, within the experimental ranges 
given in eq. (jTQJ) . For the Vl phases in eq. © we take 
fiat distributions in the interval [0, 2ir]. We make Monte 
Carlo simulations for different values of 9\ 2 and 9j 3 
mixing angles, allowing tan 2 9\ 2 an d tan 2 9 23 to vary 
respectively within their allowed intervals, in consistency 
with the lepton and quark mixing angles. From eq. ©, 
by using the fact that #13 is small and that $23 is maxi- 
mal, we get 

J ~ — sin 2#i2 sin 2#i3 sin tfi 
8 

This expression tells us that the J parameter is within 
the range |J| < 0.042. However there is a non trivial 
correlation between J and 9i 2 INS . Because the CKM 
is given by the experimental data, and V^jf is fixed to be 
zero, the phase and the 9^ NS angle are almost fixed 
as a function of 9™ NS . 

In figs. [T][3] we report the result of our simulation for 
J. We plot the correlation between the J invariant and 
sm 2 9f 2 MNS for V M BiMaximal (fig. EQ), TriBimaximal 
(fig. W), and V M with tan 2 6^ = 0.4 (fig. E). First of 
all, from fig. Q] we see that the solar mixing angle 9i 2 INS 
is constrained to have sin 2 9i 2 MNS > 0.36 for Vm Bimax- 
imal. From figs. [T][3]we see the correlation between the 
structure of Vm and the CP violating invariant J. In 
particular, for Vm BiMaximal J close to zero. For Vm 
TriBiMaximal J is around 0.042. Finally for Vm such 
that tan 2 9\ 2 = 0.4 we get that J can be any value be- 
tween —0.04 and 0.04. We also see that a better determi- 
nation of the sin 9[ 2 MNS could give a stronger prediction 
for the J invariant in the case of Vm TriBimaximal. 
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FIG. 1: The correlation between the Dirac CP violating 
parameter J and sin 2 8[ 2 MNS for Vm BiMaximal. We also 
plot the experimental central value, the la, and the 2a for 
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FIG. 2: The correlation between the Dirac CP violating pa- 
rameter J and sin 2 #f 2 MJVS for V M TriBiMaximal. 
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FIG. 3: The correlation between the Dirac CP violating pa- 
rameter J and sin 2 6>f 2 MJVS for Vm such that tan 2 6^' = 0.4. 



In figs. (US] we report the result of our simulation for 
Si (5*2 plots have similar shapes). The expressions in eqs. 
([S]) give us the range for these invariants: 

1*9x1 < 0-14 \S 2 \<0AI (14) 

We plot the correlations between the Si invariant with 
respect sin 2 9[ 2 MNS for V M BiMaximal (fig. gj), TriBi- 
Maximal (fig. 0, and such that tan 2 6\ 2 = 0.4 (fig. [§]). 
From the figures we obtain that for Vm BiMaximal the 
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FIG. 4: Same as fig [T] (Vm BiMaximal) for the correla- 
tion between the Majorana CP violating parameter Si and 
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FIG. 5: Same as fig [2] (Vm TriBiMaximal) for the correla- 
tion between the Majorana CP violating parameter Si and 
sin tf 12 
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FIG. 6: Same as figEKYu such that tan 2 8^' = 0.4) for the 
correlation between the Majorana CP violating parameter Si 
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Majorana CP invariant S\ is close to zero, for Vm TriBi- 
Maximal S\ is around 0.13. Finally for Vm such that 
tan 2 OX" = 0.4 we obtain that Si can be any value be- 
tween —0.14 and 0.14. Similar results hold for the other 
Majorana CP violating invariant S%. We see that also 
in this case a better determination of the 9i 2 MNS mixing 
angle will give a stronger constraint for the Si (and S2) 
invariant for Vm TriBiMaximal. As for J, these correla- 
tions of Si (and S2) with respect to 9 P 2 MNS are predic- 



tions of any theoretical GUT model that gives a relation 
of the type V M = U CKM Q U PMNS with V$ = 0. In 
the next section we will show how to construct an explicit 
model that predict V^ = 0. 



III. A TOY MODEL 

In this section we will show how to construct a toy 
model that gives us the relation V M = U CKM Q U PMNS 
with V^ = 0. We will not take care of explicit values 
for the masses. To obtain them we need an unwanted 
fine-tuning. However, as long as we reefer to the mix- 
ing angles only, this model can be seen as a toy model 
explaining the relationship between the CKM and the 
PMNS mixing matrix and the appearance of a zero (1, 3) 
entry in the quark-lepton correlation matrix V M . 



A. V M in theories with see-saw of type I 

Let's fix the notations in the lepton sector. Let Yi 
be the Yukawa matrices for charged leptons. It can be 
diagonalized by 



Y = C/ ; Y ; A Y/ 



(15) 



Let be Mr the Majorana mass matrix for the right neu- 
trino and Md the Dirac mass matrix. Under the assump- 
tion that the low energy neutrino masses are given by the 
see-saw of Type I we have that the light neutrino mass 
matrix is given by 



M v = M 



D 



Mr 



-M 



D ■ 



(16) 



Let us introduce Uq form the diagonalization of the Dirac 
mass matrix 



M D = UoM^vJ 



(17) 



then wc define V M by the diagonalization of the light 
neutrino mass 



M v = M 



ID 



1 

Mr 



M D = 



U v M*Ut = UqV M„(V ) Ul (18) 
Finally the lepton mixing matrix is 

jjPMNS = u \ lJv = uJ UoV M (19) 

From eqs. (|17H18[) we see that the mixing matrix V M 
diagonalize the following symmetric matrix: 



c = m£vJ-^v*m£ 



1 

Mr 



(20) 



where Vq is the mixing matrix that diagonalize on the 
right the Dirac neutrino mass matrix in eq. (|17[) . 
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B. V as correlation matrix in GUT 

In GUT models, such us generic 150(10) or Eq, there 
are some natural Yukawa unifications. These give up to 
an interesting relation between the JJ CKM quark mixing 
matrix, jj pmns lepton mixing matrix and V M obtained 
from eq. (j2"0)) . In fact V turns out to the the corre- 
lation matrix defined in eq. ([3|). In the quark sector we 
introduce Y u and Y d to be the Yukawa matrices for up 
and down sectors. They can be diagonalized by 

Y u = U u Y A Vj (21) 
Y d = U d Y d A VJ (22) 

where the Y A are diagonal and the U s and Us are unitary 
matrices. 

Then the quark mixing matrix is given by 

jjCKM = U ]jj d (23) 

In GUT models such as 5*0(10) or Eq we have intrigu- 
ing relations between the Yukawa coupling of the quark 
sector and the one of the lepton sector. For instance, 
in minimal renormalizable 5*0(10) with Higgs in the 10, 
126, and 120, we have Yi « Yj. In fact the flavor sym- 
metry implies the structure of the Yukawa matrices: the 
equivalent entries of Yi and Y d are usually of the same 
order of magnitude. We conclude that, as long as the fla- 
vor symmetry fully constraints the mixing matrices that 
diagonalize a Yukawa matrices, we have Ui ~ U d *. 
From eq. (fT9|) we get 

u PMNS ~ v^u v M 

where V is the mixing matrix that diagonalize the ma- 
trix C of eq. (|2TJ|) . If we call Y u the Yukawa coupling 
that will generate the Dirac neutrino mass matrix Md, 
we have also the relation 

Y V *Y* -> U cV: (24) 

This relation, together with the previous one, implies 

u pmns ~ v*v:v m 

If the Yukawa matrices are symmetric, for example in 
minimal renormalizable 5O(10) we have only small con- 
tributions from the antisymmetric representation 120, 
the previous relationship translates into a relation be- 
tween U PMNS , U CKM and V M . In fact we have 

Yu = Y l u -> V: = U u (25) 
Yd = Y* - VZ = U d . 
The first relation tell us that 

u PMNS = viu u v M . 

Finally, by using the second relation in eq. (|25p and the 
definition of the CKM mixing matrix of eq. (9J we get 
that 

jjPMNS _ /jjCKM\]yM 



The form of V can be obtained under some assump- 
tions about the flavor structure of the theory. This model 
will give for example a correlation V v with — 0. As 
a consequence of the from of quark- lepton correlation ma- 
trix V M there are some predictions for the model. For 
example the prediction for 9 p ^ /rNS of [HI and the correla- 
tions between CP violating phases and the mixing angle 
6»i2 of sect. El 



C. V from 5*3 flavor symmetry 

As we told V M diagonalize the matrix C of eq. (1201) . 
In this section we will show that a 53 flavor permutation 
symmetry, softly broken into 52, gives us the prediction 
of V$ = 0. 

The six generators of the 53 flavor symmetry are the 
elements of the permutation group of three objects. The 
action of 53 on the fields is to permute the family label 
of the fields. In the following we will introduce the 52 
symmetry with respect the 2nd and 3rd generations. The 
52 group is an Abelian one and swap the second family 
{a*l,(^)l, s Ll c l , n R , (u^r, s R , cr} with the third one 
{tl, (v T )L,bL,t L ,TR, {v T )R,b R ,t R }. 

Let us assume that there is an 53 flavor symmetry at 
high energy, which is softly broken into 52- In this case, 
before the 53 breaking all the Yukawa matrices have the 
following structure: 

fa b b\ 

Y = lb a b) (26) 

\b b a) 

where a and b independent. The 53 symmetry implies 
that (l/y/3, l/v3) is an eigenvector of our matrix 
in eq. (|26|) . Moreover these kind of matrices have two 
equal eigenvalues. This gives us an undeterminated mix- 
ing angle in the diagonalizing mixing matrices. 
When 53 is softly broken into 52, one get 

(a b b\ 

Y = lb c d) (27) 

\b d cj 

with c ss a and d ~ b. When 53 is broken the degen- 
eracy is removed. In general the 52 symmetry implies 
that (0,1/^,-1/^2) is an exact eigenvector of our ma- 
trix (|27|) . The fact that 53 is only softly broken into 
52 allows us to say that (1/V3, l/v3, 1/V3) is still in 
a good approximation an eigenvector of Y in eq. (|27|) . 
Then the mixing matrix that diagonalize from the right 
the Yukawa mixing matrix in eq. (|27p is given in good 
approximation by 

w -V2/V3 re 1/V3 \ 
«1/V6 wl/v 7 ^ 1/V2 (28) 
«1/V6 ^1/V3 -I/V2J 

where we did not prompt the phases. 
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Let us now investigate the V in this model. The mass 
matrix Mjj will have the general structure in eq. (f2"T)) . 
To be more defined, let us assumed that there is an extra 
softly broken Z2 symmetry under which the 1st and the 
2nd families are even, while the 3rd family is odd. This 
extra softly broken Z2 symmetry gives us a hierarchy be- 
tween the off-diagonal and the diagonal elements of Mo, 
i.e. &, d « a. c. In fact if Zi is exact both b and d are 
zero. For simplicity, we assume also a quasi-degenerate 
spectrum for the eigenvalues of the Dirac neutrino matrix 
as in fn\ . 

The Majorana right handed neutrino is of the form 

(a b b'\ 
M R = b c d (29) 
\b' d ej 

Because S3 is only softly broken into S2 we have that 
a w c ~ e, and b ~ b' w d. In this approximation the 
Mr matrix is diagonalized by a Ur of the form in eq. 
(f2"5)l . In this case we have that M v defined in eq. (fTc?)) is 
near to be S3 and S2 symmetric, then it is diagonalized 
by a mixing matrix U u near the TriBiMaximal one given 
in eq. 1281 The C matrix is diagonalized by the mixing 
matrix 

V M = V V V R (30) 

We obtain that Vm is a rotation in the (1, 2) plan, i.e. it 
contains a zero in the (1,3) entry. 

As shown in [77| , it is possible to fit the CKM and the 
PMNS mixing matrix within this model [8l| . 

IV. SUMMARY AND CONCLUSIONS 

In this paper we show the power of the quark lepton 
correlation in GUT models. In particular we investigated 



the correlation between the CP violating invariants and 
the mixing angle 12 in the lepton sector. To extract 
these informations we used a Monte Carlo approach to 
take into full account the presence of unknow unphysical 
phases in the definition of V . 

We obtain that 

0?s MNS = & 

J, Si , S 2 are correlated to sin 2 9[ 2 MNS 

with the theoretical input V^jf = only. The results 
are such that in the next future it will be possible to 
make cross check from the experimental evidences and 
discriminate the validity of this approach. 

To better clarify the importance of the quark lepton 
correlation we shown how a toy model, based on the S3 
flavor permutation symmetry and GUT, predicts the cor- 
relation matrix V M to have a zero (1,3) entry. In this 
model V M , defined by [3] is also related to the Dirac and 
Majorana neutrino mass matrix as in eq. 1201 For this 
reason its form is given by the symmetries of the model. 
In this model we have Vj^f = 0. We can apply all the 
nice predictions obtained in section UI1 
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